Enhanced electromagnetic transition dipole moments and 
radiative decays of massive Majorana neutrinos thanks to 
the seesaw-induced non-unitary effects 



Zhi-zhong Xing * and Ye-Ling Zhou t 

Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China 



Abstract 

In an extension of the standard electroweak theory where the phenomenon 
of lepton flavor mixing is described by a 3 x 3 unitary matrix V, the electric 
and magnetic dipole moments of three active neutrinos are suppressed not 
only by their tiny masses but also by the Glashow-Iliopoulos-Maiani (GIM) 
i-Ch ' mechanism. We show that it is possible to lift the GIM suppression if the 

canonical seesaw mechanism of neutrino mass generation, which allows V to 
be slightly non-unitary, is taken into account. In view of current experimental 



9* 

constraints on the non-unitarity of V, we find that the effective electromag- 
netic transition dipole moments of three light Majorana neutrinos and the 
rates of their radiative decays can be maximally enhanced by a factor of 
0(1O 2 ) and a factor of 0(1O 4 ), respectively. This novel observation reveals an 
intrinsic and presumably significant correlation between the electromagnetic 
properties of massive neutrinos and the origin of their small masses. 
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The standard model (SM) of particle physics contains three species of massless 
neutrinos which only take part in weak interactions. Since 1998, a number of convincing 
neutrino oscillation experiments have established the fact that neutrinos are actually massive 
and lepton flavors mix in a way similar to the quark flavor mixing [1] . This breakthrough has 
lots of profound and far-reaching impacts on cosmology and astrophysics, as neutrinos are 
so abundant in the Universe that their tiny masses and flavor conversions may significantly 
affect the cosmic expansion, structure formation and many other astrophysical events [2]. 

The most popular mechanism of generating finite but tiny neutrino masses beyond the 
SM is the canonical seesaw mechanism [3] , in which the small neutrino masses are attributed 
to the existence of heavy degrees of freedom such as the right-handed Majorana neutrinos. 
In this elegant picture the 3x3 lepton flavor mixing matrix V has a striking difference from 
the 3x3 quark flavor mixing matrix in the SM: it is not exactly unitary due to small mixing 
between light and heavy neutrinos as a result of the Yukawa interactions. The deviation 
of V from a unitary matrix can be appreciable, in particular when the seesaw mechanism 
works at the TeV scale so as to avoid the seesaw- induced hierarchy problem [4] and satisfy 
the testability requirement at the Large Hadron Collider [5]. A careful analysis of current 
data on the invisible width of the Z° boson, universality tests of electroweak interactions, 
rare leptonic decays and neutrino oscillations indicates that the unitarity of V is good at 
the percent level and possible non-unitary effects are also allowed at the same level [6]. 

Now that three active neutrinos acquire their respective masses in the seesaw mechanism, 
they should have the electromagnetic dipole moments (EMDMs) through quantum loops. 
The fact that the Majorana neutrinos are their own antiparticles implies that they can only 
have the transition EMDMs between two different neutrino mass eigenstates in an electric or 
magnetic field [2]. The relevant radiative decays of the heavier active neutrinos, which may 
contribute to the cosmic infrared background in the Universe [7], are of particular interest in 
cosmology. Since all the previous calculations of the EMDMs and radiative decays of massive 
neutrinos were done by assuming V to be exactly unitary, we find it highly necessary and 
important to recalculate the same quantities by taking into account the seesaw-induced 
non-unitary effects. 

The aim of this paper is just to examine the seesaw- induced non- unitary effects on the 
EMDMs and radiative decays of active Majorana neutrinos. We do a complete one-loop cal- 
culation and demonstrate that such effects are likely to be significantly larger than the stan- 
dard (unitary) contributions, because the suppression induced by the Glashow-Iliopoulos- 
Maiani (GIM) mechanism [8] in the latter case can now be lifted. A careful numerical 
analysis shows that the effective EMDMs of three neutrinos and the rates of their radiative 
decays can be maximally enhanced by a factor of O(10 2 ) and a factor of O(10 4 ). Such an 
intrinsic and presumably important correlation between the electromagnetic properties of 
massive neutrinos and the origin of their small masses must be taken seriously, and it may 
even serve as a sensitive touch-stone for the highly-regarded seesaw mechanism in the future. 



The canonical seesaw mechanism is based on a simple extension of the SM in which 
three heavy right-handed neutrinos are added and the lepton number is violated by their 
Majorana mass term [3]: 

-C v = I^Y„HN R + \MM r N r + h.c. , (1) 
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where H = ia 2 H* with H being the SM Higgs doublet, £ L denotes the left-handed lepton 
doublet, N R stands for the column vector of three right-handed neutrinos, and M R is a 
symmetric Majorana mass matrix. After spontaneous SU(2) L ® ^(1)y ~~ ^ ^(l) e m § au g e 
symmetry breaking, H achieves its vacuum expectation value (H) = v/y/2 with v ~ 246 
GeV. Then the Yukawa-interaction term in C u yields the Dirac mass matrix M D = Y v v/y/2, 
but the Majorana mass term in C v keeps unchanged because right-handed neutrinos are 
the SU(2) L singlet and thus they are not subject to the electroweak symmetry breaking. 
The overall neutrino mass matrix turns out to be a symmetric 6x6 matrix and can be 
diagonalized through 

M M>--ro i)> <*> 

in which we have defined M v = Diagjm^ m 2 , m 3 } and M N = Diag{M 1? M 2 , M 3 } with 
and M i being the physical masses of three light neutrinos v i and three heavy neutrinos N i 
(for i = 1, 2, 3). The 6x6 unitary matrix U is decomposed as [9] 



U J\S BJ\0 

where 1 denotes the 3x3 identity matrix, U and V are the 3x3 unitary matrices, and A, B, 
R and S are the 3x3 matrices which characterize the correlation between the active or light 
neutrino sector (V ) and the sterile or heavy neutrino sector (U ). A full parametrization of U 
in terms of 15 mixing angles and 15 CP- violating phases has been given in Ref. [9]. One may 
express the flavor eigenstates of three active neutrinos in terms of the mass eigenstates u i and 
iVj. In the mass basis of charged leptons and neutrinos, the leptonic weak charged-current 
(cc) and neutral-current (nc) interactions read 



(V m u lL + R ai N iL ) W~ + h.c] , 

2 cos 9„ 



"£ nc = {W^V)^ + N ih r(rfRhN jh + [^(V^R) tJ N 3L + h.c] } , (4) 



where a runs over e, /i or r, V — AV Q is responsible for the flavor mixing of active neutrinos 
u i , and R measures the strength of charged-current interactions of heavy neutrinos (for 
% — 1, 2, 3) [10]. A small deviation of V from V is actually characterized by nonvanishing R, 
as VV"t = AA^ = 1 — RR) holds. The exact seesaw relation between the masses of light and 
heavy neutrinos is VM U V T + RM N R T = 0, which signifies the correlation between neutrino 
masses and flavor mixing parameters. 

Let us consider the radiative u i — > + 7 transition, whose electromagnetic vertex can 
be written as 

r&(o) = M«(i ^0 + ^(^9,75) (5) 

for a real photon satisfying the on-shell conditions q 2 = and = 0. In Eq. (5) and 
/i^ are the electric and magnetic transition dipole moments of Majorana neutrinos, and their 
sizes can be calculated via the proper vertex diagrams in FIG. 1 (weak cc interactions). The 
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^-Z self-energy diagrams in FIG. 2 (weak nc interactions) do not have any net contribution 
to e { j and but we find that they play a very crucial role in eliminating the infinities 
because the divergent terms originating from FIG. 1 are unable to automatically cancel out 
in the presence of the seesaw- induced non- unitary effects (i.e., R ^ and V ^ V Q ) unless 
those divergent terms originating from FIG. 2 are also taken into account. This observation 
is new. It implies that the non-unitary case under discussion is somewhat different from the 
unitary case (i.e., R = and V^V = VqV = 1) discussed before in the literature [11], where 
the Feynman diagrams in FIG. 2 are forbidden and the divergent terms arising from FIG. 1 
can automatically cancel out. 

After a careful treatment of the infinities and non-unitary effects in our calculations [12], 
we arrive at 



ieG T 
eG p 

— = — I m { — m.- 



(6) 



where 



Fa 4 



2-6, _ ?L gglnj. 

i-C (i-o 2 a-O' 



(7) 



with £ a = m^/Myy (for a = e, /x, r) denotes the one-loop function. Although this result is 
formally the same as that obtained in Ref. [11], they are intrinsically different as the seesaw- 
induced non-unitary effects on (j,^ and were not considered in the previous works. To see 
how important such non-unitary effects may be, let us make two analytical approximations. 
First, F a ~ 3 (2 — £ Q ) /4 holds to a good degree of accuracy for £ a <C 1. Second, V = AV ~ 
V — TV is also a good approximation for small non-unitary corrections to V^,, where [9] 
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(8) 



with c,„- = cos and s,-,- 



iS- 

e « 



sin#j- (here 0^ and 5^- are the mixing angles and CP- 
violating phases). Note that the light-heavy neutrino mixing angles 8 ik (for % — 1,2,3 and 
k — 4, 5, 6) are at most of (9(0.1) [6], such that the deviation of V from is at the percent 
level or much smaller. Then we obtain 



E F a {V^) ~ -| £ ( T K))aj + (TV )«i (^] - 7 E [*« (KL (K)^ 



(9) 
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The first and second terms on the right-hand side of this equation correspond to the non- 
unitary and unitary contributions, respectively. While the former is suppressed by sj k ^$ 
(9(1CT 2 ) (for i — 1, 2, 3 and k — 4,5, 6) hidden in T, the latter is suppressed by £ a < 4.9 x 1CT 4 
(for a = e, ji, r) due to the GIM mechanism. We therefore draw a generic conclusion that 
the seesaw-induced non-unitary effects on and /i^ can be comparable with or even larger 
than the standard (unitary) contributions. 

The rates of radiative v i — >■ i/ - + 7 decays are more sensitive to the non-unitarity of V, 
since they directly depend on |/x^-| 2 and | e^- 1 2 . Namely, 



with /i cfr = y |/iy| 2 + I e^- 1 2 for v i — > +7 being the effective EMDMs and /i B = e/ (2m e ) be- 
ing the Bohr magneton. The size of V , can be experimentally constrained by observing 
no emission of the photons from solar v e and reactor V e fluxes. More stringent constraints 
on /i eff come from the Supernova 1987A limit on neutrino decays and from the cosmological 
limit on distortions of the cosmic microwave background radiation (in particular, its infrared 
part): /i cff < a few x 1CT 11 /i B [13]. Now that more and more interest is being paid to the 
cosmic infrared background relevant to the radiative decays of massive neutrinos [7], it is 
desirable to evaluate /j, s and V , , on a well-defined theoretical ground, such as the 

i j 

canonical seesaw mechanism under discussion. 



3 We proceed to numerically illustrate the non-unitary effects on /x ff and T u , , . 

i 3 

In view of current neutrino oscillation data [14], we take Am^ ~ 7.6 x 10~ 5 eV 2 , Am 2 2 ~ 
±2.5 x 10~ 3 eV 2 , 9 12 ^ 34°, 9 23 ~ 45° and 6 13 ~ 9° [15] as our typical inputs. The mass 
scale of three neutrinos, the sign of Am 2 2 and the values of three CP-violating phases of V 
remain unknown. In our numerical calculation we consider both normal (Am 2 2 > 0) and 
inverted (Am| 2 < 0) hierarchies by fixing the smallest neutrino mass to be 5 meV, and allow 
all the CP-violating phases to vary between zero and 2n. Those small active-sterile neutrino 
mixing angles in Eq. (8) are constrained by present experimental data [6] as follows: 

T u < 5.5 x 1(T 3 , |T 21 | < 7.0 x 10" 5 , 
T 22 < 5.0 x 1(T 3 , |T 31 | < 1.6 x 10~ 2 , 

T 33 < 5.0 x 10~ 3 , |T 32 | < 1.0 x 10" 2 . (11) 

Note that all s ik in T (for % = 1,2,3 and k = 4,5,6) are positive or vanishing. The CP- 
violating phases 8 ik are all allowed to vary from zero to 2n, but they must satisfy the above 
constraints together with the relations VV^ + RB) = 1 and VM U V T + RM N R T = 0. It is 
in principle possible to determine the masses of three heavy Majorana neutrinos with the 
help of the exact seesaw relation, if all the other parameters are known [16]. Although such 
a prediction is in practice impossible, we find that the upper bound on the non-unitarity of 
V implies the existence of a lower bound on M v In other words, the values of M i cannot be 
too small (i.e., they should be far away from those of so as to avoid too significant non- 
unitary effects. To assure that radiative corrections to the masses of three light neutrinos (via 
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the one-loop self-energy diagrams involving the heavy neutrinos) are sufficiently small (e.g., 
smaller than 0.5 meV) and stable, we simply assume that the masses of three heavy neutrinos 
are nearly degenerate [17] and not more than 0(1) TeV. This assumption implies that we are 
concentrating on a limited and safe parameter space of the TeV seesaw mechanism, but it 
is instructive enough to reveal the salient features of the non-unitary effects on the effective 
EMDMs /i e ff( z/ i ~> v j + 7) and the radiative decay rates r 

To present our numerical results in a convenient way, let us define 



p — 



X! ( S lfc + S lk + S lkj 
.fc=4 



1/2 

(12) 



which measures the overall strength of the unitarity violation of V, and e uv e [0,0.15) is 
reasonably taken in our calculations. Namely, we allow each s ik (for % — 1, 2, 3 and k — 4, 5, 6) 
to vary in the range < s ik < 0.15. The numerical dependence of /i eff (^ — > Vj + 7) and 
r . , on e uv is shown in FIGs. 3 and 4, respectively. Some discussions are in order. 

(1) Switching off the non-unitary effects (i.e., e uv = 0), we obtain the effective electro- 
magnetic dipole moments 
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(14) 



for the inverted mass hierarchy with m 3 ~ 5 meV, where the uncertainties mainly come 
from the unknown CP- violating phases 5 12 , S 13 and 5 23 - Such standard (unitary) results are 
far below the observational upper bound on /i eS (< a few xlO -11 ji B [13]), but they serve 
as a good reference to the non-unitary effects on fi eS being explored in this work. 

(2) FIGs. 3 and 4 clearly show that /i eS and r v _^ v +y can be maximally enhanced by a 
factor of C(10 2 ) and a factor of (9(10 4 ), respectively, in particular when e uv approaches its 
upper limit as set by current experimental data. Note that the magnitude of A* e fj(z / 2 ~~ ^1+7) 
may be strongly suppressed in the inverted neutrino mass hierarchy. The reason is rather 
simple: on the one hand, m l ~ m 2 holds in this case, and thus e 12 oc (m 2 — m x ) must be 
very small; on the other hand, fj, 12 depends on lm(V al V* 2 ), so it can also be very small when 
the CP- violating phases are around zero or ir. This two- fold suppression becomes severer 
for the decay rate r +7 , because it is proportional to (m 2 — m 1 ) 3 / u 2 ff (z/ 2 — > u 1 + 7). 

(3) Our numerical analysis shows that the results of /i s and T , are sensitive to the 
absolute neutrino mass scale for both normal and inverted mass hierarchies. For instance, 
A*eff ( u 2 ~ * v \ + 7) an d fJ-esi^ ~ * v \ + 7) § e t enhanced when m l changes from zero to 5 meV 
in the normal mass hierarchy; while /x eff (i / 1 — > v 3 + 7) and /i cff (^ 2 — > v 3 + 7) are enhanced 
when m 3 changes from zero to 5 meV in the inverted mass hierarchy. This kind of sensitivity 



6 



is not so obvious if one only takes a look at the expressions of /i^ and in Eq. (6). The 
main reason is that a change of m l or m 3 requires some fine-tuning of the active-sterile 
neutrino mixing angles and CP-violating phases as dictated by the exact seesaw relation 
VM U V T + RM N R T = 0, leading to a possibly significant change of // eff . The dependence of 
T v + on the absolute neutrino mass scale is somewhat more complicated, as one can see 
fromEq. (10). 

(4) We find that the CP-violating phases play a very important role in fitting both the 
exact seesaw relation and Eq. (11). If the heavy neutrino masses M i are not suppressed, 
then an appreciable value of e uv requires some fine cancellations in the matrix product 
RM N R T such that sufficiently small m i can be obtained from VM V V T = —RM N R T . On 
the other hand, we remark that it is actually unnecessary to require M i to be around or 
above the electroweak scale. The seesaw-induced non-unitary effects on /i eff and r . , 

i j 

can be significant even if one allows one, two or three heavy neutrinos to be relatively light 
(e.g., at the keV mass scale). Such sterile neutrinos are interesting in particle physics and 
cosmology. Note that it is easier to satisfy the exact seesaw relation with an appreciable 
value of e uv by arranging M i to lie in the keV, MeV or GeV range. This kind of low-scale 
seesaw scenarios [18] might be technically natural, but they have more or less lost the seesaw 
spirit. Of course, sufficiently large M i and sufficiently small 6 ik can always coexist to make 
the seesaw mechanism work in a natural way, but in this traditional case the non-unitary 
effects are too small to have any measurable consequences at low energies. 

It is also worth pointing out that the seesaw-induced non-unitary effects on ■ and 
are rather different from the case of making a naive assumption of the flavor mixing between 
three active neutrinos and a few light sterile neutrinos [19]. The latter can directly break 
the unitarity of the 3x3 active neutrino mixing matrix V and then lift the GIM suppression 
associated with /i^ and e^. This kind of non-unitary effects are not constrained by the 
seesaw relation, and thus they are more arbitrary and less motivated from the point of view 
of model building. 



4 1 We have explored the seesaw-induced non-unitary effects on the electromagnetic 
transition dipole moments and radiative decays of three active neutrinos. We find that such 
effects are possible to be comparable with or larger than the standard (unitary) contribu- 
tions, because the suppression induced by the GIM mechanism in the latter case can be 
lifted. Our numerical analysis has illustrated that the effective electromagnetic dipole mo- 
ments of three neutrinos and the rates of their radiative decays can be maximally enhanced 
by a factor of O(10 2 ) and a factor of (9(10 4 ), respectively, no matter whether the seesaw 
scale is around or below the TeV energy scale. This observation is new and nontrivial, and 
it clearly reveals an intrinsic and presumably important correlation between the electromag- 
netic properties of massive neutrinos and the origin of their small masses. Such a correlation 
may even serve as a sensitive touch-stone for the highly-regarded seesaw mechanism. 

We thank W. Chao, Y.F. Li and S. Zhou for useful discussions. This work was supported 
in part by the National Natural Science Foundation of China under Grant No. 11135009. 
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FIGURES 




FIG. 1. The one-loop Feynman diagrams (and their charge-conjugate counterparts) contribut- 
ing to the EMDMs of the Majorana neutrinos, where a = e,fi, r and i, j = 1, 2, 3. 




FIG. 2. The one-loop j-Z self-energy diagrams (and their charge-conjugate counterparts) as- 
sociated with the EMDMs of massive Majorana neutrinos, where / denotes all the SM fermions 
and i,j = 1, 2, 3. 
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FIG. 3. Illustration of the seesaw-induced non-unitary effects on n eS for three active neutrinos. 
The standard (unitary) results correspond to e uv = 0, and their uncertainties come from the three 
unknown CP-violating phases of V . 
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FIG. 4. Illustration of the seesaw-induced non-unitary effects on T 
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neutrinos. The standard (unitary) results correspond to £ uv = 0, and their uncertainties come 
from the three unknown CP-violating phases of V^. 
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